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Abstract 



We consider various processes of single top production at LHC in the 
theoretical framework of the MSSM and examine the role of the supersym- 
metric electroweak one-loop corrections in a special moderately light SUSY 
scenario, in an initial parton-pair cm. high energy range where a logarith- 
mic asymptotic expansion of Sudakov type can be used. We show that the 
electroweak virtual effects are systematically large, definitely beyond the 
relative ten percent size, particularly for a final tH~ pair where a special 
enhancement is present. We show then in a qualitative way the kind of 
precision tests of the model that would be obtainable from accurate mea- 
surements of the energy distributions of the various cross sections and of 
the top polarization asymmetries. 
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I. INTRODUCTION 



The process of single top production in proton-proton collisions has recently been 
described in detail [1] and its relevance in theoretical models of electroweak physics has 
been stressed with special emphasis. In the particular framework of the Standard Model, 
the relevant final pairs can be (at partonic level): 1) a top and a light quark (t-channel 
production), 2) a top and a b (s-channel production), 3) a top and a W~ (associated 
production), and the corresponding Feynman diagrams at Born level are shown in Fig. 1,2. 
As stressed in Ref. [1], from precise measurements of these processes the first direct 
determination of the CKM Wtb coupling at hadron colliders will be obtained, which 
already motivates the existence of dedicated experimental and theoretical studies. 

The theoretical analysis of the first two processes has been performed in a number of 
papers, both within the SM and beyond it, with special emphasis on the production at 
TEVATRON [2]. The third process will only be measured at LHC, and to our knowledge 
a detailed theoretical study at the perturbative one-loop level, valid in particular for the 
MSSM, does not exist yet. 

In this paper, we shall assume a preliminary discovery of supersymmetry at LHC 
(perhaps already at TEVATRON) and a SUSY scenario of a "moderately" light kind, 
in which all sparticle masses lie below, say, approximately three-four hundred GeV. Our 
aim will be that of showing that, under these conditions, single top production at LHC 
might provide crucial, accurate precision tests of the candidate theoretical model. For this 
preliminary analysis we shall stick to the MSSM, but our treatment could be modified in 
a straightforward way to examine less simple supersymmetric proposals. 

In our analysis, we shall consider "large" values of the initial parton pair cm. energy 
y/s, typically of the 1 TeV size. There exist two main reasons that motivate our choice. 
The first one is that within this energy range we shall feel entitled to make use, at the 
one-loop perturbative level, of a simple asymptotic logarithmic expansion of so called 
Sudakov kind, whose validity in the MSSM is related to the assumption of a light SUSY 
scenario, so that s/M 2 (the typical parameter of the logarithmic expansion) is of order 
ten (M is the heaviest SUSY mass of the scenario). The second one is that, for y/s values 
much larger than the final t, W masses, remarkable kinematical simplifications arise whose 
theoretical origin will be shown. These make, in particular, the treatment of the associated 
production process much simpler, and allow to understand the potentialities of its relevant 
supersymmetric effect from inspection of short analytic formulae. 

The plan of this paper will be the following. In Section 2A we shall briefly derive the 
expressions of the electroweak Sudakov expansions at one loop to next-to leading order, i.e. 
retaining the squared and the linear logarithmic terms, for the t- and s-channel processes. 
For the latter, a partial MSSM calculation of the Yukawa effect is already available, in 
an arbitrary cm. energy configuration where in general all the parameters of the model 
will enter the theoretical predictions [3]; our calculation will include logarithmic terms of 
several origins, but owing to our energy choice a strongly reduced number of parameters 
will remain in the relevant expressions, as we shall show in detail. Section 2B will be 
devoted to a more accurate description of the determination of the Sudakov expansion for 



2 



the so called associated tW production process. Here, as we said, a theoretical calculation 
in the MSSM is not known to us; independently of this, we shall comment anyway a few 
special features of our expansions, in particular the simplifications due to the choice of 
a high energy regime. Within the MSSM, another process to be added to the previously 
considered ones seems to us to be that of tH~ production. This has also been already 
studied in an arbitrary energy configuration [4], [5], [6]; we shall derive here in Section 
2C only its electroweak logarithmic Sudakov expansion, to be compared with the three 
previous ones to have a more complete description of the asymptotic single top production 
at LHC. 

A final topic that seems to us to deserve some consideration is that of the single 
top polarization. This is automatically fixed in the t, s channel processes in the chosen 
high energy region where the final top must be of left-handed type; this feature is not 
valid for tW~ and tH~ production. As a consequence, one may in principle consider 
longitudinal polarization asymmetries for these processes. In Section 2D, we shall derive 
the asymptotic expressions of these quantities, that will turn out to be relatively simple. 

Having at disposal all the relevant asymptotic expressions for the four considered 
processes, we shall then devote Section 3 to an investigation of which kind of precision 
test of the model could be provided by realistic experimental measurements. This will be 
done in a necessarily qualitative way, exploiting the available preliminary conclusions of 
Ref. [1], assuming that the distributions in y/s can be determined with a certain accuracy 
in our chosen range. Our presentation will remain for the moment at a qualitative level. 
In fact, the purpose of our paper is mostly that of showing the possible relevance of 
measurements that appear to us, at least in principle, performable, and to encourage an 
experimental performance like that indicated in Ref. [1]. This point will be discussed in 
the final conclusions, given in the short Section 4. 

II. SINGLE TOP PRODUCTION 

As we said, this Section will be devoted to a derivation of the asymptotic electroweak 
expressions for the four considered processes. To save space and time, since the major 
part of the relevant technical definitions and properties has been already exhaustively 
illustrated in previous papers [7], we shall quickly write the final formulae for the (t,d), 
(t,b) and (t,H~) processes, that have already been studied in the literature (but not in 
our chosen scenario) as we anticipated. The treatment of the final (t, W~) pair will be 
slightly more detailed since for this process, as we said, we are not aware of a theoretical 
investigation in the MSSM, and also since in this case the benefits of the chosen high 
energy configuration deserve, in our opinion, a short mention. 
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A. T,S CHANNELS SINGLE TOP PRODUCTION 



1. bu — > td with W exchange in the t channel 

We begin with the t-channel process represented in Fig. la at Born level. As one 
sees, there is only one diagram with W exchange in the t-channel, t — (p b — p t ) 2 . The 
corresponding scattering amplitude is: 

AB ° rn = ^TT + ^[u(d,T')rPMu, \')}[u(t,r)^P L u(b, A)] (2.1) 

where X,X',t,t' are the b,u,t,d helicities, Pr l = (1 ±7 5 )/2 are the projectors on R, L 
chiralities ans s w is the sine of the Weinberg angle. 

It is convenient to work with helicity amplitudes F\ t \i jT y; retaining only the top mass 
and setting all the remaining masses equal to zero leaves one single amplitude F : 



.Born 47iaSy/P 



(2.2) 



with (3 = £ = 1 - ^. 

The expression of the differential cross section is: 



dcosO 11527rs 

spinfiol 

which becomes after color average: 

da Bom f3 2 na 2 s 



<h ; E \Fxyml 2 (2.3) 



dcosO 8s^/(t — myy) 2 



(2.4) 



At one-loop, the Sudakov electroweak corrections can be of universal and of angular 
dependent kind, and we follow the definitions of [7]. 

The effect of the universal terms on the helicity amplitude can be summarized as 
follows: 

F Univ_ = pBorn_ ]_ j + c ™{uu) L + C eW [dd) L + C eW (tt) L \ (2.5) 

where: 

c ew (qq)L = c ew (qq) L = c(qq, gauge) L + c(qq, yuk) L (2.6) 

- ail + 26cL) s o s 

c{dd, gauge) L = c{uu, gauge) L = —— — 2 2 (2 log — log — r ) (2.7) 

L4:4:7rs w c w fi^w m w 
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c(bb, yuk) L = c(tt, yuk) L = -— r [log—] [-f(l + cot 2 (3) + -f(l + tan 2 f3)\ 

(2.8) 

Note that we use m w as a scale for the log term. This is mandatory for the squared 
log which arises from W loops. For what concerns the single log, as we do not consider 
"constant terms" this a matter of choice and we decide to keep m w as a scale for the 
above terms. 

The angular dependent terms have the following expression: 

pang _ pBorn j ^png j (2 9) 



c-. = -^^±^[.o g ^l[lo g ^] 



a(l - 10c 2 w ) -t 

2 2 log jl- o 



2 2 [log— ] [log -j-] (2.10) 



At high energy we have t ~ — f (1 — cos#) and u ~ — f (1 + cos#). 

There are also SUSY QCD corrections, only of universal kind. We shall treat them 
as "known" terms for what concerns our search of electroweak effects; at one loop, their 
expression is not difficult to be derived, and reads: 

pUniv SUSYQCD = pBorn j og S j ( 2 11) 

3tt M| usy 

In addition to the previous terms of Sudakov type, there are at one-loop "known" 
linear logarithms of RG origin, whose expression we quote for completeness: 



1 r a ^dF 



Born 



F-_ = -^^p][lo gw 



w 

a 2 sy/]3 



using the lowest order Renormalization Group (3 function for the gauge coupling g = 
e/s w : ft* = - \ in MSSM. 



2. ud — > tb with W exchange in the s channel 



Neglecting again all quark masses with the exception of m t one finds from the Born 
diagram of Fig. lb: 



I Born 



2S 2 V {S - 77%) 



[v(d, \')rPMu, A)] [u(t, r) la P L v{b, t')\ (2.13) 
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There are now in principle two helicity amplitudes F\,a',t,t' : 

(2.14) 



,Bom _ 2nasy/J3(l + COS0) 



s 2 w {s-m 2 w ) 



^ = _ (2 . 15) 

s 2 w {s-m 2 w ) 



r 

where 3 = 1 — ^ 

s 

Note that the second amplitude F^°™ + is depressed with respect to the first one 
by the factor m t /y/s, so that in our chosen configuration it will be treated with some 
approximations, neglecting its one loop corrections. 

In the differential cross section the two contributions will be: 

1 a = To^ X! |-^A,A',r,r'| 2 (2-16) 

da Bom (3 2 na 2 s r/1 m2 m 2 . 2 „. ,^ _ 

^ = 32^(.-ma,)» ' (l + CMg) ' + T°"'«l < 2 ' 17) 

where the first term comes from F_ + h and the second one comes from the "depressed" 

amplitude that will be retained in Born approximation. For the relevant amplitude 

we obtain the following effects: 

l 
2 

with the expressions of the various c ew given for the t-channel case and 

= a2{ \ + / OS6) [\o g -^]{[A(Q dQb + QuQt) + ^ 2 +f^ ] log ^] 

-[KQdQt + + gig%t|^g4 ] log Z^ ]} (2.19) 

s w c w s 

with <?f L = 2/3 - 2s^Q,. 
We list also the SUSY QCD logarithms: 



F Univ + = F Born + _ [ + C ^( MM ) L + C e -(dJ) L + C e »(tt) L ] (2.18) 



- log m|^; ] (2 - 20) 



and those of RG origin: 



1 ~ AF Born q 
F R U = ~ T^^iS^] log- ' 



4vr 2Ly " rf^ 2 J b M 2 



w 



2sUs-m 2 w ) l ° g M 2 w {2J1) 
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using (3° = — \ in MSSM and choosing M w as the RG scale, which seems to us 
justified at our logarithmic level. 

For the t- and s-channel processes, the logarithmic effects at one-loop are fully sum- 
marized by eqs. (2.5-2. 12, 2.18-2.21). They represent an original result of this paper, and 
provide an effective simple representation in the scenario that we have chosen , valid, we 
repeat, to next-to leading logarithmic order. As one sees, the welcome feature of the 
formulae is that in the logarithmic expansion there appears only one SUSY parameter, 
tan (5. This would allow a remarkably simplified test of the model, that will be discussed 
in the next Section 3. 



B. ASSOCIATED t,W PRODUCTION: b g ->■ t W~ 

We now move to the process to which we devote, for various reasons, a more detailed 
description in our scheme, i.e. the so called associated (t, W) production. Fig. 2 shows 
that now there are two Born diagrams, one (a) in the s-channel with exchange of a bottom 
quark and one (b) in the u = (pb — Pw) 2 channel with exchange of a top quark. Denoting 
by p = \pt\ = \pw\ the modulus of the final cm. momentum, we define: 



2 P o'o 2 Pw P 



(yfs-m t ) 2 - mfy 



y/s H y/s n E t + m t ^ (y/s + m t ) 2 - 
and the W coupling -gwh^pPh^v) w hh g WL = 



m 2 w 



(2.22) 



Neglecting all the quark masses with the exception of the top one, one gets the helicity 
amplitudes F x ^ x >,n> with the helicities A(6) = ±1/2, fi(g) = ±1, X'(t) = ±1/2, fi'(W) = 
p! = ±1 for transverse or p'(W) = for longitudinal Wq respectively. 



jpBama _ 1 n „ ( X \ . / Wj + m tY ~ m W r r 

\,n,\',n' ~ 2 9wl9s ^)\I ~ ° X < L i 

1 cos - [1 + pp cos 6 + 2A(/i + p cos 9) + 2p'(p + 2A) sin 2 -]5 x ,y 

Y t 

— -sin-[2A(l + pp'cos6) + p + p! cos 6 - 2p'(l + 2Xp) cos 2 -]5x,-\> } 

(2.23) 



F Borna_ n f A' V(>/i + Wj) 2 - < 

{ ^p(p + 2X)sm e -(P + p'°)S x , x , 

+±^±(l + 2\p)cos 9 -((3-p'°)6 x ,_ x , } (2.24) 
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pBornb 9WL9sS , \\ ( y/s + m t ) 2 - m 2 w 



- p'°)(l + /i/i'(2 cos 9-1)- 2A(/i + //)) 
-(1 - cos0)(/x// + - ////(/? + cos6>) - (1 + (3 cos 9) 
+2A(/i(l + (3cos9) + n'tf + cos 9) - (1 - cos 0)(// + cos -<S A>A , 

H -{(1 - P' )(-2/j,' cos 2 - + // + // cos 6> - 2A(2 / u/i / cos 2 - + 1 + /i/i' cos 0)) 

2 2 2 


-/i(l + /3cos0) -//(/3 + cos0) -2cos 2 -(/i/3 + /i') -4A(/i// + /3) cos 2 - 

+2A////(/3 + cos 0) + 2A(1 + /3 cos 0)} sin -5 A ,-v 

+ i/i^ x ' x '^r 1 cos 2 (_1 + - 2 cos ^ + 2A(/U + ^ 

1 + r 6 1 

+5A,-A'^- J sin-(/i' + 2A(1 + /i/i' + 2////cos0))] } (2.25) 



A' S\ (Jl + m t ) 2 — 
F w , - 9wl9s(j)( 4Mw{u _ ml)V z 



L 



{ sin ^(1 + cos 0)2/3[/i(/3 - 0° + 1) + 2A]5 A ,y 

cos -2/3[(l + 2A/i cos 0)(1 + (3 - + 2A/i + cos 6]6 x ,-\> 

+ 7^ A ' A '^ sin 2^'° (2A " M1 + 2 cos ^ " P(2X + ^ 

+5A,-v^^cos^(-/?'°(l-2A/i(l-2cos0)-/3(l + 2A/i))] } (2.26) 

where A' are the color matrices associated with the gluon and g s is the QCD coupling 
constant. 

With the color sum J2i ^ r [(y)(y)] = 4 the cross section is 



2^2 

da (3 



dcos9 7687rs eo 

spins 



E I^W^I 2 (2.27) 



In our special scenario we are allowed to neglect m t / 'y/s, m w /y/s, but not m t /M w . 
This leads to a remarkable simplification of the previous expressions, that become now: 



F%3J - 9wl9s{^) cos ^(1 + 2A/i) V (2.28) 



F?°™J - 9wL9skj) cos ^(2A/iCOS0 - 1)* A>L (2.29) 
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- sin 2 (/i + 2A) (2 - 30) 



^Sv.o - "^l^(^)5a,l{ -j^- sin % + 2\)8 XtX , + cos + 2Ap)tf A ,_A' 



V2m t l + 2A/icos0 

-M^^- A,COS 2 [ l + cos0 ]} (2 ' 31) 

Explicitly, the only remaining amplitudes are 

F Born_a +b ^ (2 32) 

2 cos | 

F Bo™ +0+ 6 ^ ^ s (^) 2 cos | (2.33) 

„ A' V2m f 1 - COS 

f ++0 ^^(y)i^r cos 2 ( TT^0 ) (2 - 34) 

They provide a rather simple expression of the differential cross section: 

d ° B ° rn _» _ ma - fj 2 , u 2 , ZM^I f2 35) 

rfcos0 24a^t« 2la + + 2M^ J ^ 

The electroweak Sudakov terms at the one-loop level are now relatively simple to 
compute and to show. More precisely, we obtain for the universal component for transverse 
W T production: 

= F^LJ \ ( c™(bb) L + c™(ti) L ) + c™{W T )] (2.36) 

where c(bb), c(tt) have been already defined by eqs. (2.6-2.8) (following the notations 
of [7,8]) and 

c-(^H£H-log^] (2.37) 

and for longitudinal Wq production: 

<t + ,o = F-™,ol \ ( <^(»)r + ^(^R ) + c ^(^o)] (2.38) 
where in MSSM 
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(2.39) 



such that 

s ir 13 + 14c^ 

ji 

For the electroweak angular terms we find: 



= [ - s ]|log 4 l{ [log vit** + St' + i log v > 

^- B r-,J-^][>o g ^]{[.o^i[^f] + ^.o g ^} 

(2.41) 



- *f£W- ][log w}{ [- log -] - — ^ log -T } (2.42) 



a iri s irr 4 -£ 1 - lOcL , — u 
— H lo g77^]{ b log—] "a" 

•'W 1V1 W ° 6 5 T4/ 

There are also SUSY QCD universal terms: 

jpUniv SUSYQCD _ jpBorn r lri _ g 1 ( A o\ 



SUSY 



jpUmv SUSYQCD = ^ } _^_] (3.44) 

~'+' + '° -.+,+.01 3?r S M S 2 USY J V ; 

For this process, there are no one loop ew RG terms. 

A few comments are now appropriate. The first one is that we have obtained a great 
lot of simplification neglecting systematically terms of order m t / \fs in the components of 
the helicity amplitudes. This is in agreement with the general strategy of our asymptotic 
scenario, that assumes \fs to be sufficiently larger than all the masses of the process 
and will retain only the two leading logarithmic terms of the Sudakov expansion. For a 
different scenario, e.g. one where is smaller, the remaining components of the cross 
section must be retained and the theoretical formulae will be more complicated. Their 
computation might be necessary if for instance the (optimistic?) assumption of light 
SUSY were not verified. Then the region of large ^/s would have no special reasons (of 
Sudakov origin) to be priviledged, and thus a complete one-loop calculation would be 
imperative, to be used for instance in a relative low y/s range where a Sudakov expansion 
would certainly not be acceptable. This can certainly be done, but in our optimistic 
attitude we will now show which interesting features might arise from a measurement of 
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the process at LHC, if our chosen scenario were ( kindly) chosen by Nature. The second 
point concerns Eqs. 2. 34, 2. 38, 2. 39, 2. 42. As one can verify, they are in full agreement with 
the equivalence theorem which states that, neglecting M^/s terms, the expressions for 
longitudinal Wq production are identical to those for the Goldstone boson G~. This 
can be seen by taking the expressions written in the next Section 2C for tH~ production 
and replacing Mfj with Mq = M$y, m t cot/3 with m t , and rattan/? with -m^, and this 
constitutes a positive check of our calculation. 

In fact, having concluded the study of the production in the MSSM of a top and a W 
boson, it seems almost natural to enlarge the treatment so as to include a process which 
is tightly connected with it, i.e. that of production of a top and of the charged Higgs of 
the model. This process has been already studied [4-6]. In particular, in [5] one can find 
an exhaustive discussion of the SUSY QCD corrections, particularly relevant in the large 
tan/? limit. In the next Subsection 2C we shall show that the purely electroweak super- 
symmetric corrections considered in [4], but not treated in [5], are also very important 
especially in the electroweak Sudakov regime of the LHC in which we are interested. 



C. ASSOCIATED t, H~ PRODUCTION: bg -► tH~ 

This process is similar to the longitudinal part of the previous process bg — > tW~ . 

There are two Born diagrams like in Fig.2, but replacing the final W~ by the H , (a) 
s-channel with bottom quark exchange and (b) w-channel with top quark exchange. We 
use the same notations as before with u = (p b — Ph) 2 and define: 



n = tyt o>o = ^h_ = Pt = (V^-m t y-m 2 H 
1 V~s P ~ V~s 1 E t + m t \(^- s + mt )i-ml 

and the H~tb coupling ~ e [mt cot (3Pl + mj, tan (3Pr\ . 

Neglecting the quark masses, except m t , one gets the helicity amplitudes Fx^y with 
the A(fe), //(<?), X'(t) helicities respectively: 



.Born a _ _^_ ( g / (v^W ~ ™jr { 2X( l^l ) 

w s w M w y 2 J \ s x y 2 ' 2 



+( 2 - )cos-8x-\> }5x^[m t cot/35x,L + m b tan/35x, R } (2.46) 



~* s w M w (u-m?) { 2l s { ^^[(^-)sm-2AV 

+ (^—) cos -5x,-y\ 

sin^ A , A ,[-(2A)p t (l + cos 9)5 X ,^ + (2A)(p t cos0 - E t + V7s)5 x J 

I _|_ f t Q 

H-^ 1 ) cos ^x,-y[Pt(l - cos 9)5 x ,-n + (Pt cos 9- E t + Vs)S x J } 

[m t cot f35 x ,L + m b tan f35 x , R ] (2.47) 
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With the color sum X^(y)(y) = 4 the cross section is 



da 



dcos9 768ns 



(2.48) 



spins 



neglecting m t / \/s,m H / 'y/s, but not m t /m w nor m b /M w we have 



rpBom a+b 



eg s \/s X 1 9 l — cos9 n<r 

[ — )cos-(— )dx-yOx-^[m t cot Pdx >L + m b t&n(3d\, R \ 



s w M w 2 2 v l + cos6>' 



(2.49) 



explicitly 



rpBorn a+b 
*-,+,+ 



eg s m t cot (3 . \\ 9 1 — cos 9, 

( — ) cos — ( 

s w M w y 2 J 2 v l + cos#' 



(2.50) 



rpBorn a+b 



eg s m,b tan (5 . A' 9 \ — cos9 
s w M w ^~2> C ° S 2^1 + cosr 



(2.51) 



da Born naa s (m 2 cot 2 /3 + m\ tan 2 /3)t 2 

dcos6> 48s 2 y M^ / ws 2 



(2.52) 



Moving to one loop, one finds in our scenario: 
a) One loop electroweak universal terms 



f~:+:+ = f-%:+\ - ( ^ ew m L + <r^) R > + c ™ +t+ {H-)\ 



(2.53) 



?Univ 



F*°?L[ l - ( c ew (bb) R + c ew (ti) L ) + <%_ t _(H-)] 



(2.54) 



with 

:+ah- 



l + 2c 2 



-fh 



^log 2 



+IW _?_][- 17 +10c ^ + 



3m 2 

( 1 +tan 2 /5) + - x -^(l + cot 2 /5)]} 



4s 2 y M 2 , 



(2.55) 



-+ 



a; 



1 + 2c^y 2 



8s 2 r 2 



log^ 



5 + 22c 2 



+ [1 ° g M^ ][ 3Qs 2 w c 2 w l.s 2 V 2 



+ 



m * (l + cot 2 /3) + -^ r (l + tan 2 /5)]} 



■>W 1V± W 



4s 2 yM 2 / 



(2.56) 
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such that 



^^H[-.o g ^]||±i|k]} (2.57) 



£]{[-. og^][^|f]} (2.58) 
b) One loop electroweak angular terms 



n:i- - f??a- nog - s^t' > (2 - 60) 



Again, there are no one loop ew RG terms, 
c) One loop SUSY QCD universal terms 



a. 



-pUniv SUSYQCD _ -pBorn r s , - i /<-, fi -,x 



^ StVSyQCD = pBorn r_ «. j (2 g 2) 

3tt 8 M| usy J 



As shown in [9] for the decays t — > iJ + — > t6 and in [4,5] for this bg — > 

process, in addition to the SUSY electroweak and QCD contributions coming from loop 
diagrams, there are important renormalization terms in the large tan f3 limit which affect 
the Born value Eq.(2.50), namely one should replace mbtan/3 by TO;,tan/3/(l + A&). The 
explicit expression of A& in terms of the various mass and mixing parameters is given 
in [5] . This correction becomes important when masses are not degenerate and tan (3 is 
large, so that we will keep track of this term in the applications that we shall make in 
Section 3. 

The (long) list of equations that we have derived in this three subsections provide the 
relevant asymptotic expressions of all the cross sections of the four considered processes 
of single top production. They represent an original result of our paper, that cannot be 
compared with other one loop formulae since these do not exist, to our knowledge. Note 
that the Standard Model result can be easily reproduced by our equations with the simple 
formal rules that have been already discussed in previous references and that correspond 
essentially to the following replacements: for the gauge part: 21og(s/M^) — log 2 (s/Mf^) 
-> 31og(s/M^) - log 2 (s/M^), and for the Yukawa part: 2m 2 (l + cot 2 /?) -> m\ and 
2m 2 (l + tan 2 (3) — > m 2 in the heavy quark case; m 2 (l + cot 2 /5) — > m 2 and m 2 (l+tan 2 (3) — > 
ml in the Higgs and Goldstone cases. 
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In our study we have not considered until now the possibility of defining different 
observables, e.g. exploiting the fact that the polarization of the final top should be, in 
principle, observable [1]. This might lead to the definition of longitudinal polarization 
asymmetries already considered in the case of top - antitop production [10]. As a matter 
of fact, such observables should not provide exciting information (apart from possible 
tests of the V-A structure) for the first two t, s-channel processes, where at high energy 
the final top is supposed to be essentially of left handed type. But for the two remaining 
processes this property is a priori no longer valid, and the production of a right handed 
top is not forbidden. This is made possible by the nature of the other particles produced 
in association (i.e. either a longitudinal W boson or a scalar charged Higgs), that breaks 
chirality conservation. We have thus decided to devote the next Section 2D to a collection 
of the relevant formulae for the related polarization asymmetries. Again, we remark that, 
to our knowledge, these formulae do not exist in the literature. Also, a preliminary dis- 
cussion of the expected experimental and theoretical uncertainties of their measurements, 
which exists for the four unpolarized production cross sections [1,11] is not yet available. 
Therefore our formulae should essentially be considered as a theoretical proposal. The 
reason why we decided to show them is that we believe that they exhibit some interesting 
features, that would deserve, in our opinion, a further investigation performed in a more 
realistic approach. 



D. top polarization asymmetries in t, W and t,H production 

We define in general 

_ 4a l lo °v{bg -> t L Y~) da 1 lo °v{bg -> t R Y~) 

at[e) ~ 1 d^d I^e J/ 

M l ^(bg^t L Y-) da 1 iQ°P(bg^t R Y-) 

[ dcosO dcos6 J 1 ' ' 

where Y = W or H . 

In the two considered processes we shall have 
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(2.65) 



The asymptotic expansions at one loop of the various helicity amplitudes that enter 
equations (2.61-2.62) have been already derived in Sections 2B-C and therefore we shall 
not write the one loop expressions of the asymmetries. It might be useful, though, to 
observe that a t {6) in both cases is not vanishing at Born level. More precisely, we obtain 
the following expressions 



14 



s 2 + u 2 



tW~ : 4 Born) (9) = — — 2 -^— (2.66) 



S 2 + U 2 + ™ 

W 



tH~ : 4 B ° rn \0) = i tan ^"i CQt ^ (2.67) 
v ' m\ tan 2 (3 + m 2 cot 2 f3 v ' 

Note that both asymmetries are energy independent at the Born level. Their behavior as 
functions of cos# (for final tW~) or tan/5 (for final tH~) are shown in Fig. (3). 

We have now concluded the list of the asymptotic Sudakov expansions of the considered 
observables of single top production. All the formulae have been given at the partonic 
level. In the following Section 3, we shall try to derive more realistic predictions for the 
observable physical processes. 



III. APPLICATIONS TO OBSERVABLE PROCESSES 



A. Unpolarized cross sections 



We shall begin our analysis with the investigation of the electroweak one-loop effects 
in the four considered unpolarized cross sections. With this aim, we shall first provide 
a calculation of the inclusive differential cross section of the actual processes, defined as 
usual as (Y = d, b, W~ or H~): 

da(PP^tY + ...) 1 f™se ma da^ty 

ds = S Le min dC0s9 [ ^ ^ (T ' C0S ^^W (S) ] (3 ' 1} 

where t = f , and (ij) represent the initial partons of each process. are the corre- 
sponding luminosities 

L y (T,oose)= r aX dy [i( x )j(l)+j( x )i(L)] (3.2) 

where S is the total pp cm. energy, and i(x) the distributions of the parton % inside 
the proton with a momentum fraction, x = e 5 , related to the rapidity y of the tY 
system [12]. The parton distribution functions are the 2003 NNLO MRST set available 
on [13]. The limits of integrations for y can be written 

1 1 

y max = max{0,min{r - -log*, Y + -logx, -log(^/r)}} 

ymin ymax 

(3.3) 

where the maximal rapidity is Y = 2, the quantity x is related to the scattering angle in 
the tY cm. 
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1 -|- cos 6 



(3.4) 



X = 



1 — cos 9 



and 



cos 6. 




s 



(3.5) 



expressed in terms of the chosen value for pT,min which gives the integration limits for 
cos 6 in eq.(3.2). 

The results of our analysis are shown in Fig. (4) (energy dependence at various tan/3) 
and (5) (tan/3 dependence at two special energies y/s = 0.7 and 1 TeV). The value of 
PT,min is always chosen to be 50 GeV. As one sees clearly from the Figures, there are two 
main general features that emerge, i.e.: 

a) the one-loop SUSY electroweak effects are systematically large (well beyond the 10 
% level) in the considered energy range, particularly for the t channel process and for the 
tH~ process, where they can reach the 40% level for tan/3 = 50. In this last case, they 
are as important as the SUSY QCD renormalization effects considered in [5]. 

b) The effects are strongly dependent on tan/3. They have a minimum value for 
tan/3 ~ 10, then they increase regularly for either smaller or larger values. 

Encouraged by this preliminary information, we have tried to perform a y 2 fit to 
tan /3 assuming a reasonably realistic experimental determination of the \fs distribution, 
which also takes into account expected theoretical uncertainties. With this purpose, we 
have proceeded in the following way. First, we have separated the first three processes of 
single top production (i.e td, tb and tW~) from the fourth tH~ case. For the first three 
cases we have followed the pragmatic attitude of assuming, from the general conclusions 
given in [1], that the measurements of all the cross sections can be performed in the 
energy range 500 — 1500 GeV (with 20 GeV binning) with an overall (theoretical and 
experimental) accuracy of 10 %. In fact, an ambitious final goal of 5% was mentioned in 
these conclusions, so that we might say to have followed a reasonably conservative attitude. 
With this overall uncertainty, we have obtained the results of combined conventional x 2 
fit that are shown in Fig. (6). In more details, we define for Y = d,b, W~ 



as in Eq.(3.1). Then, for each true value tan/3* of the unknown tan/3, we compute 

2 / + R , ^ ^ ^y(g,tan/3)-0 y (g,tan/3*) V 

X tan /3, tan /3 =2^ 7 7 ^ > 3 - 7 

t^V cry(s, tan/3*) J 

where the sum over s runs over the above range of y/s and the assumed experimental 
error a Y is fixed by the request of 10 % accuracy. Of course, x 2 (tan/3*, tan/3*) = 0. 
We now vary tan/3 until we have x 2 (tan/3, tan/3*) = 1. This determines two values 
tan/3 = tan/3* + A± with A + > and A_ < shown in Fig. (6). 



CV(s,tan/3) 



da 
ds 



(PP^tY +•••), 



(3.6) 



16 



As one sees from Fig. (6), the determination of tan/3 performed in this way might be 
remarkably successful, in particular for the large tan (3 ^ 20 region, which is expected to 
be poorly determined at the LHC time [14]. Here the error of the fit will be reduced, 
under our working assumptions, to the small 2-3% value for tan (3 ~ 50. Although, we 
repeat, our preliminary analysis has been undoubtedly qualitative, we believe that this 
result could be considered as a first encouraging step. 

Before presenting our results for the (t,H~) production process, a few preliminary 
remarks seem to us to be appropriate. The first one is that, using the notations of Ref. 
[5], we have limited our investigation to the "inclusive" process, i.e. gb — ► tH~ (the same 
we did in this paper when we considered the inclusive tW~ production). In Ref. [5] this 
process has been determined to NLO in QCD, including the SUSY QCD corrections, for 
a general cm. energy configuration. The results of the analysis show that the dominant 
higher order SUSY QCD correction is not due to virtual loops, but rather to an effect 
of "coupling renormalization" type, that can be described by a constant (cm. energy 
independent) parameter defined as A b (eq.(7) of that reference), that depends on several 
MSSM parameters. For certain choices of the latter (like e.g. large gluino masses) the 
effect of At, on the cross section can be rather large; more precisely, it can produce 
modifications in the 30 — 40 percent range that cannot evidently be neglected and must 
be accurately taken into account in any realistic fit to the data. 

The second remark that we want to make is that the calculation of the virtual elec- 
troweak SUSY corrections to the inclusive (tH~) process at high energy that we have 
performed is not appearing in Ref. [4,5]. We repeat again that our analysis is truncated 
at the next-to leading logarithmic order, where only the first two (quadratic and linear) 
terms of a logarithmic expansion are retained, since for a preliminary indicative investi- 
gation we believe that this approximation should be sufficient. 

After these remarks, we are now ready to show the results of our investigation. Briefly, 
we assumed an expression for the inclusive cross section, written in analogy with eq.(3.1), 
in which the SUSY corrections contain both the dependent component of Ref. [4,5] 
and our logarithmic corrections of Sudakov origin. The correction A& actually depends on 
several MSSM parameters. An inspection of Table I in [5] shows that the largest effects 
are obtained in a region of the parameter space that roughly corresponds to large values 
of tan (3. In such an extreme case, a simple parametrization of is 

A 6 = 0.004 tan (3. (3.8) 

The coefficient 0.004 is fixed in order to generate a negative correction to the cross section 
of approximately (relative) forty percent for tan/3 = 50 as shown in Table I of [5]. Of 
course, our choice is rather drastic, but we used it to show the qualitative features of the 
competition between this term and our electroweak components. Adding the electroweak 
Sudakov correction, we have repeated the \ 2 analysis performed in the previous case of 
the three (t, d), (t,b), (t,W~) processes. 

More precisely we have first simulated actual data for da/ds in the process bg — > tH~ 
by using our theoretical expression with the full set of logarithmic Sudakov corrections as 
well as the term A5. These fake experimental data are assumed to have an experimental 
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error at the fixed value of 20 %, a reasonable choice as supported in [5]. Again, we 
considered values of y/s between 0.5 and 1.5 TeV with 20 GeV spacing. As a second step, 
we made a x 2 determination of tan j3. With this aim, we have made a fit to the simulated 
experimental data with different theoretical representations labeled "A b ", "e.w.", and 
"e.w.+A fe " according to the set of corrections that are taken into account. The meaning 
of the three choices is the following: 

1. da th /ds = da Ab /ds. We examine what would be the bias and the confidence bounds 
in the determination of tan/3 when the electroweak Sudakov logarithms are not 
included in the theoretical expression. 

2. da th /ds = da e ' w '/ds. As before, but retaining in the theoretical expression for the 
cross section only the electroweak Sudakov logarithms. This support the statement 
that the correction A& is necessary in any realistic investigation in points of the 
parameter space where the extreme parametrization (3.8) applies. 

3. da th /ds = da e - w - +Ah jds. In this case, we use as a theoretical representation, the 
expression with which we generated the simulated experimental data. Of course, 
there is no bias in this case and the best estimate of tan/3 coincides with the true 
value used to simulate the experimental data. Here, we just compute the limits and 
compare them to those obtained in the above two cases. 

Hence, following the notation of Eq. (3.6-3.7), we have considered, for this analysis, 

O h (s, tan (3) = ^(PP^t #" + •••), (3.9) 
ds 

»- lw r-£( ^y w )' (-) 



s 



where "kind" can be "A 6 ", "e.w.", and "e.w.+A fe " as discussed above. Due to bias in the 
cases "kind" = "e.w." and "A 6 ", the minimum of x 2 as a function of tan/3 at fixed tan/3* 
is not reached at tan /3 = tan /3* and also it is not zero, Xmm > 0- As before, we vary tan (3 
until Ax 2 = X 2 ~Xmm = 1- This determines a confidence interval tan /3~ < tan /3 < tan /3 + . 
In the left part of Fig. (7), we show the curves for tan^ for the various choices of the 
theoretical representation of the cross section. In the right part of the same Figure, we plot 
instead the percentual value of the ratios (tan P ± — tan/3*)/ tan/3* in order to emphasize 
the bias. Of course, in the case "kind" = "e.w. + A b " , there is no bias at all. 

From the analysis of Fig. (7), several comments arise. The various stripes obtained 
with different theoretical representations of the cross section have approximately the same 
width, i.e. a ±1—2% with respect to the central value. The bias that is obtained when one 
of the two kinds of radiative corrections is not included in the theoretical representation 
of the cross section is about —10% for tan/3 = 15 and increases up to about —20% 
at tan/3 = 50. These values match the corresponding typical values of the corrections 
themselves in the cross section. 
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It seems that we can conclude that the two kinds of corrections, in this extreme 
case, have similar effects and neglecting one of them would lead to a large error in the 
estimated central value for tan /3, even if the confidence interval remains small. Of course, 
neglecting at all the two corrections would be definitely unacceptable. If, on the other 
hand, we accept that the sum of corrections "e.w. + gives a reliable representation 
of the data, then we see that it is possible to estimate tan (3 with a rather small error in 
the range 2-4 % for all the explored values of tan/5 between 15 and 50. This remarkable 
feature is clearly a consequence of the fact that the cross section depends on tan (3 already 
at the Born level. 

Of course, one must remember that the precise form of the correction A& depends on 
several MSSM parameters that are in principle not known. In practice, with no a priori 
knowledge on the parameters entering A& (i.e. sbottom squark masses, gluino mass, /i 
parameter, . . .), the above determination of tan/3 is expected to give a tight confidence 
interval of a few percent with a systematic bias which is roughly of the same order of the 
partially or totally unknown correction due to A&. 



B. Polarization asymmetries 



To conclude our analysis we present now the results of a determination of the 1- 
loop SUSY Sudakov effects on the two polarization asymmetries defined at parton level 
by eqs. (2.61,2.62). For this preliminary determination following the treatment for the 
unpolarized case we defined the inclusive asymmetries 

*M = + ^ <*■"> 

with 

d<J L ,R 1 [ c ° s 0ma* nJVi^t^nY 



ds S 



/ dcos9 [ £ L^r, cos 9) r L f {s) ] (3.12) 

J COS u m i n ■■ a cos u 

i j 

Figs. (8,9) show the effect of the one loop corrections at variable energy and variable 
tan p. From these figures one can derive the following indications: 

a) in the (t,W~) case the effects varies, changing sign with tan/5, moving from a 
positive value of 2 percent for tan (5 — 2 and reaching a negative value of approximately 
5 percent for tan/3 = 50. 

b) in the (t, H~) case, the situation is rather peculiar. The main feature that emerges 
is that the value of the asymmetry at Born level is already tan (3 dependent and reaches 
the limiting values ±1 in correspondence to small (^ 6) or large tan/5 values. In this 
ranges, one obviously finds that the one loop effect is practically vanishing. But the 
interesting fact that one observes is that, looking at the sign ±1 of the asymmetry, one 
would be able to fix the proper range of tan f3 with extreme reliability. This property of 
the inclusive (t,H~) polarization asymmetry has not been, to our knowledge, stressed in 
the literature and it seems to us that it would deserve a more detailed investigation. 
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IV. CONCLUDING REMARKS 



The basic assumption of this paper has been a direct discovery of supersymmetry 
at LHC (perhaps already at Tevatron), in a moderately light SUSY scenario, where all 
sparticle masses satisfy the bound < M, with M of the order of ~ 350 — 400 
GeV. Under these conditions, we have considered, in the theoretical framework of the 
MSSM, four processes of single top production in the kinematical configuration y/s (initial 
partons center of mass energy) ~ 1 TeV, that will be accessible at LHC. The purpose 
of our study was that of investigating whether measurements of distributions of cross 
sections or of other observables, performed under realistic theoretical and experimental 
precisions, might lead to stringent consistency tests of the electroweak sector of the model, 
at its perturbative one- loop level limit. Our main tool has been the use of a logarithmic 
Sudakov expansions for the considered sector, which appears theoretically justified in the 
considered scenario at the chosen cm. energy. The latter expansion has been computed 
to next to leading order, i.e. retaining the quadratic and linear logarithmic terms and 
neglecting further contributions, essentially constant (energy independent) ones. In such 
an approach, the only supersymmetric parameter that enters the electroweak expansion 
is tan/3. Under these conditions, one might hope that a suitable fit to realistic data can 
fix with a reasonable accuracy the correct tan (3 value. The outcome might be, depending 
on the case, either a confirmation of a previous alternative determination (e.g. for small 
tan/3 < 20 values) or possibly a prediction for large tan (3 values (tan/3 > 20), where an 
alternative determination might still be lacking at the LHC running time. If the second 
more appealing and ambitious programme turned out to be the correct one, the outcome 
of "predicting" the tan/3 values from precision measurements at LHC would assume a 
theoretical relevance comparable, to a certain extent, to that belonging to the memorable 
prediction of the top mass from precision measurements at LEP1. 

To perform a realistic investigation, i.e. one that might lead to a measurable prediction 
for a physical process, a first necessary condition seemed to us that of examining whether 
appreciable effects would exist in the basic partonic processes treated under reasonable 
qualitative overall precision assumptions. The results of this qualitative investigations 
show that the overall logarithmic electroweak effects would be systematically large in the 
considered yfs region with relative values between 30-40 % in all four processes for large 
tan/3. Starting from this encouraging discovery we have shown that from a combined 
fit to the distributions of the cross sections of the first three processes a remarkable 
"determination" of tan/3 would be possible. For the fourth (tH~) case, we have shown 
that the situation would be less simple since a previously determined SUSY QCD constant 
correction to the Born amplitude has to be considered, that would depend on several other 
parameters of the model. This would make the fit less simple, although still potentially 
crucial for the large tan/3 values, for which we believe that this type of searches would 
be particularly relevant. In this respect, we have also investigated, although in a more 
qualitative an essentially preliminary way, the possible extra information obtainable from 
measurements of top polarization asymmetries, finding results that, a priori, would seem 
encouraging, particularly for the previously discussed (tH~) inclusive production process. 
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To make the final move to realistic measurable predictions requires now a few final 
steps that we list and that we considered, though, beyond the purposes of this first 
preliminary investigation. The first step would be to transform our theoretical one-loop 
predictions, given (as it is normally done) in terms of the initial partons cm. energy y/s 
into predictions given in terms of the final tY pair invariant mass M t y. This requires 
investigations of the features of, e.g., the possible gluon emission from the final state 
that are, in principle, performable. In fact, at the moment an analysis of this kind for 
the case final top antitop production is being carried on [10]. The second step would be 
that of computing realistic theoretical and experimental uncertainties for the considered 
inclusive distributions. In [1,5] a preliminary analysis was performed for all four processes 
but a rigorous updated systematic and complete determination of the theoretical and 
experimental uncertainties does not exist to our knowledge. 

A third step would be the extension of our Sudakov expansion to (at least) the next 
to next to leading order. This requires the determination of a possible constant (e.g. 
energy independent) term which might depend on several parameters of the model. An 
example of this statement has been provided by the consideration of the constant At 
term in the (tH~) process which depends on the sbottom squarks and gluino masses, on 
/i and A b . Note though, that this A b term is of strong SUSY origin. In our case, we thus 
expect a priori that a possible constant electroweak SUSY term should be depressed by, 
roughly, a ~ a/a s factor with respect to A b which explains why we feel that it should not 
affect appreciably our preliminary results, that are due to squared and linear logarithmic 
enhancements. 

As a matter of fact, and as we anticipated at the beginning of this paper, we consider 
our long list of equations and our presentation of qualitative figures not really as final 
precise predictions, but rather as a proposal for realistic near future detailed studies. In 
this sense, we hope that this paper might soon be considered as the starting point of 
future successful and relevant supersymmetry studies at LHC. 
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FIG. 3. Born values of the final top quark polarization asymmetry for the two processes 
bg — > tW~ or tH~ . Notice that in the first case the asymmetry is independent on tan/3, but is 
scattering angle dependent. In the second process, these features are reversed. 
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FIG. 5. Percent correction to the (unpolarized) cross section do Ids due to the full set of 
Sudakov logarithmic corrections. We show the curves in the MSSM as functions of tan (3 at two 
reference energy values yfs = 0.7, 1 TeV. 
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FIG. 6. Results from the x 2 analysis of the tan (3 dependence of the cross sections for the pro- 
cesses bu — > td, ud — > tb and bg — ► tW~ . A 10% precision is assumed on all the three processes. 
For each hypothetical true value of tan /3, the figures show the Ax 2 = 1 boundaries tan [3 + A± 
(shown in the inset plot) and the corresponding relative error defined as ^(A + — A—)/ tan/3. 
The Figure is cut at the upper limit of 25% since for larger values of the accuracy, the present 
derivation of the confidence bounds deserves a more detailed treatment. Actually the accuracy 
is minimum around tan/3 ~ 6 and is again better than 25% when tan/3 < 3. To avoid confusion 
we have shown the results of the analysis for tan (3 > 20. 
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FIG. 7. Results from the \ 2 analysis of the tan (5 dependence of the cross sections for the 
process bg — > tH~ . We show the confidence limits tan/? 1 * 1 (defined in Sec. IIIA) on tan/3 and 
the corresponding relative bounds (tan (3 ± — tan (3) / tan (3 for different choices of the theoretical 
representation of the cross section as discussed in the text. In more details, the left figures plots 
the two curves corresponding 
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FIG. 8. Energy dependence of the final top polarization asymmetry. In the top line, we show 
for the two processes bg — ► tW~ , bg — ► ii7 _ the Born values of the asymmetry. In the bottom 
line, we show the percentual correction due to the full set of Sudakov logarithmic corrections. 
We show the curves for the Standard Model and for the MSSM at three reference values of tan (3. 
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FIG. 9. tan/3 dependence of the final top polarization asymmetry. In the top line, we show 
for the two processes bg — ► tW~ , bg — ► ii7 _ the Born values of the asymmetry. In the bottom 
line, we show the percentual correction due to the full set of Sudakov logarithmic corrections. 
We show curves for the MSSM at two reference energy values ^fs = 0.7, 1 TeV. 
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